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Abstract: A detailed study of the photophysical properties of several 
novel color centers in chemical vapor deposition diamond is presented. 
These emitters show narrow luminescence Unes in the near infra-red. 
Single photon emission was verified with continuous and pulsed excitation 
with emission rates at saturation in the MHz regime, whilst direct lifetime 
measurements reveal excited state lifetimes ranging from 1-14 ns. In 
addition, a number of quantum emitters demonstrate two level behavior 
with no bunching present in the second order correlation function. An 
improved method of evaluating the quantum efficiency through the direct 
measurement of the collection efficiency from two level emitters is pre- 
sented and discussed. 
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1. Introduction 

The emerging field of quantum optics has established a demand for an accessible solid state 
system which can generate a stream of single photons on demand 111. Although single pho- 
ton emission has been demonstrated from quantum dots (QD) HI 121, single molecules |[3| and 
nanowires|j4]|, the operation of those systems is often limited by the temperature. Diamond crys- 
tals, on the other hand, offer the most promising platform for generation of robust, photo stable, 
single photons at room temperature |I5]|6|. Nevertheless, out of more than 500 existing optical 
centers in diamond |21, only three known centers demonstrate single photon emission, namely, 
the nitrogen-vacancy complex fSjISI, the silicon-vacancy complex ||9l and the nickel-nitrogen 

complex Idol iniiia. 

Recent progress in materials science and fabrication techniques of diamond color centers 
ifTSl unveiled novel quantum emitters operating in the MHz regime lfT4l[T5l . These new centers, 
are comparable to QD in terms of brightness 1 1 1, but have the tremendous advantage of stable 
operation at room temperature. These break throughts have opened a new avenue to investigate 
essentially an unknown group of single photon emitters in diamond. 

The simple two-level nature of a select few single photon emitters allows a direct estimate 
of the collection efficiency associated with the optical setup, the precise knowledge of the col- 
lection efficiency can lead to more accurate estimations of the quantum efficiency of the more 
common three level single photon emitting systems. This type of source can thus eventually 
be used as a reference for an absolute characterization of the quantum efficiency along with 
more precise estimations of the collection efficiency of the system, provided the photo-physics 
is well understood. One can envisage that single photon sources could be implemented as a 
"single photon standard", able to link classical radiometric measurements to the fundamental 
quantum optical entities. In particular in the long term such a single photon standard could 
contribute to a re-definition of the standard units for optical radiation in terms of the "quantum 
candela" 1161 . The controlled fabrication of the MHz class of diamond single photon emitters, 
together with a recent emission enhancement by coupling the light to plasmonic structures ifTTl 
or cavities [TSl, is expected to progress the technology beyond the point in which they can 
become practical quantum information and metrology devices. 

In this paper, a detailed study of the photophysical properties of the novel family of sin- 
gle photon emitters reported in ifTsl has been undertaken. To gain insight into the electronic 



level structure of the emitters, their fundamental optical properties such as lifetime, photo- 
luminescence (PL) spectrum, sub-poissionian statistics and saturation count rate have been 
studied and compared. We concentrate on the typical PL lines in the range of 740-770 nm, 
which exhibit single photon emission with count rates ranging from L3-3.2 10^ counts/s at 
saturation. 



2. Theoretical background 

The physics of the photoluminescence emitted from a single diamond color center can be de- 
scribed within the framework of a two level system comprising a ground and excited state, and 
in many instances as a three-level system ifSl fTTlfTgl . whereby the single center is excited from 
its ground state to the excited state, with a third longer lived state providing an additional decay 
path from the excited state, as seen in Fig.[T] The rate equations describing the populations of 
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Fig. 1. Schematic diagram of a (a) two-level and (b) three-level system, where 
1,2,3 are the transition rates from level {n) to level (m). 



the two and three level systems can be written in matrix form as: 
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By solving Eqs ([T]) and (|2| with the initial condition /ii = l,n2 = 0,n3 = (i.e the system is 
prepared in the ground state), the instantaneous emission probability of a photon being emitted 
from the excited state, n2(f), can be obtained. By normalizing «2(f) to the probability of a 
photon being emitted at an infinite time «2(°°) an analytical expression describing the second 
order correlation function can be obtained. The respective analytical expressions for the two 
and three level llTTlfTgl cases are then given by: 

g(2)(T) = i_exp(-AiT). (3) 

where Xi = r 12 + I'll and 

g(2)(T) = 1 - (1 +fl) exp(-AiT) +fl exp(-A2T), (4) 
where A2 = + r23'"i2/Ai, and a = ri2r23/(Air3i). 

The decay rate Xi depends on the excitation optical power, Popt as: 

Ai =rfi(l + aPopt), (5) 



where is the inverse of the excited state Ufetime (T21) and a is a fitting parameter. 

Similarly, for the three-level system, the decay rate Ao, assuming constant with the 
excitation power ifTTI . can be written in terms of r'^^ in the limit of zero optical power 

where j3 is a fitting parameter. 

With an accurate determination of the second order correlation function one can determine 
the individual decay rates involved in the system and gain an understanding of the energy level 
structure. Upon inspection of the two and three level auto-correlation functions there is a clear 
distinction in the shape and nature of the exponential component. In the two level case the 
function exhibits a simple exponential rise at a rate equivalent to the fluorescence decay rate 
of the excited state in the limit of zero optical excitation i.e. ri2 0, with the exponential 
asymptotic to a g^^'(T) value of 1. The three level expression, on the other hand, contains two 
exponential components each with a characteristic time constant. Depending on the transition 
rates from the excited state to the shelving state and shelving state back to the ground state, 
g(^'(T) can increase beyond 1 for times > T21, before decaying back to 1 at times >> T21. 

The phenomenon of g'^'(T)^l is commonly termed "photon bunching" and this behavior 
enables a clear distinction to be made between two and three level systems. The second order 
correlation function describes the probability of detecting a photon with a delay time T after one 
photon has been detected at T=0. Hence, the bunching effect describes the enhanced probability 
to detect a photon at short times then at longer times. Indeed after the system undergoes a 
transition to its shelving state, which has longer lifetime, there is a longer interval between 
the photons. Once the system relaxed to the ground state and undergo a full emission cycle 
again, the normal photon rate is achieved again. The waiting interval while the system in the 
metastable state thus creates the bunching effect of the g'^'(T). 

The measurement g'^^\T) is done typically using a Hanbury-Brown and Twiss (HBT) in- 
terferometer whereby the coincidence events of two detectors are measured. In practice, the 
temporal jitter of the detectors and of the electronics can be ignored provided that the life- 
time is much longer than the temporal jitter; with the identification of emitting centers with 

increasing short fluorescence lifetimes Igl fTTlfTSjl . the effect of the temporal jitter on the meas- 

(2) 

ured auto-correlation function gmeas(f) must be taken into account. Therefore the measured 



auto-correlation function is given by the convolution of Eqs ( 1314b with the instrumental time 
response function J{t): 

gSL(T)- r dT'g(^\x')J{x-T'). (7) 



The majority of single photon emitters identified in diamond (such as: Si-V |[9j, NV l[20| and the 
NES lTOllTll ) have been attributed to three level systems which suffer from quenching due to 
the presence of a shelving state. However, the recent identification of highly efficient two level 
systems in diamond |fT3l [Ml has attracted considerable interest, particularly due to its central 
role in determining more accurately the collection efficiency of the setup and the quantum 
efficiency of three level systems. The quantum efficiency, tJqe, of an emitter is defined as the 
probability of an absorbed pump photon resulting in emitted photon. For a two level emitter, 
the quantum efficiency is equal to 1 when ri2 >> ^"21, indicating that for each absorbed pump 
photon a photon is emitted from the excited state. 

We can deduce the number of emitting photons or the fluorescence count rate of a two level 



system by solving the rate equations in Eq. 
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where (j) represents the single photon count rate, <j)oo = '"21'?' is the saturation count rate for 
^opt °° where 77 represents the total collection efficiency and Psat = 0!^' is the optical satura- 
tion power 

In the three level case the quantum efficiency can be derived analytically in terms of the 
detected fluorescence rate by solving the rate equations in Eq. ©I 



3. Experiment 

The CVD diamond crystals employed in this work were grown to an average size of few hun- 
dreds nanometers from diamond seeds (4-6 nm) on a sapphire substrate using a microwave 
plasma enhanced CVD technique (900 W, 150 Torr) li2l1 . A home built confocal microscope 
with a spatial resolution 400 nm and HBT interferometer were used to identify the emitting 
centers and measure the time correlation of photoluminescence (PL) intensity (Fig. 2(a)). A 
fiber coupled CW diode laser emitting at 682 nm and a 690 nm pulsed diode 200 ps pulse width 
(repetition rate from 10 MHz up to 80 MHz) were interchanged for excitation. The lasers polar- 
ization was controlled by a Glan Taylor polarizer and a half-wave plate. The diamond sample 
was mounted on a piezo XYZ stage with 0.2 nm resolution, allowing 100x100 jxin^ scans. The 
unwanted residual laser line was eliminated by a dichroic beam splitter and a Fi broad (10 nm 
FWHM) band-pass filter centred at 740 nm, 749 nm or 760 nm depending on the PL from a 
specific crystal. The PL from the emitting centers was then coupled into a 62.5 /im core multi- 
mode fiber, which acts as an aperture. A 50:50 fibre coupled beam spUtter guided the photons 
to two single photon counting detectors (APDs) and their outputs were sent to the start and stop 
inputs of the time correlator card. The PL was recorded without the bandpass Fi filter using a 
fiber-coupled spectrometer with a cooled CCD array. All the measurements were performed at 
room temperature. 

4. Results 

More than 10 different sapphire substrates were scanned, with crystals grown to a size of a few 
hundred nanometers. Fig. 2b shows a typical confocal map of the diamond crystals obtained 
by 682 nm laser excitation. The PL spectra of the bright spots revealed emission with zero 
phonon lines (ZPLs) centered at 744±2 nm (FWHM ~ llnm),749±2 nm (FWHM - 4nm), 
756±2 nm (FWHM - 1 1 nm), and 764±2 nm (FWHM ~ 10 nm), as shown in Fig. 3. In some 
cases two or more of these lines were found in one crystal, as shown in the inset of Fig. 3. 
Anti-bunching measurements on crystals containing multiple emission lines generally showed 
multiple emitters with auto-correlation functions at the zero delay time of ^'^•'(0) > 0.5, with 
the intensity of the lines varying from crystal to crystal. These characteristic emission lines 
are in agreement with cathodo-luminescence lines from chromium related centers in diamond, 
exhibiting narrow lines in the region of 740-770 nm | l22l . 

Figure 4(a-d) show the antibunching behavior of the PL lines centered at 744 nm, 749 nm, 
756 nm and 764 nm measured with the HBT interferometer for excitation powers below and 
above the optical saturation power P^at- The dip at zero delay time indicates a single photon 
emitter The raw coincidence data were corrected for the background as described in ref. ||6ll . 
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Fig. 2. (a)Experimental setup. A diode laser at 682 nm operating in CW and a pulsed 690 
nm laser were alternatively coupled to a single mode fiber, and then collimated (C) and 
polarized by a Glan Taylor polarizer (GTP). The laser polarization was varied by an half- 
wave plate (X /2). A variable neutral density filter (NDF) was used to change the excitation 
power Samples (S) were excited by focusing the laser light by a high numerical aperture 
(0.9) objective with lOOx magnification. A dichroic mirror (DM)transmitting from 720 nm, 
was used to separate the laser line form the sample fluorescence emission, when collected 
back from the same Obj. F12 are band-pass filters, 794±80 nm and alternatively 760±12 
nm or 740±12 nm, to isolate the single photon emission lines. A 100 mm focal length 
lens was used to send the single photon emission to a multimode-fiber, providing an aper- 
ture for the confocal imaging. Finally a 50:50 fiber beam splitter was used to verify the 
single photon emission by performing the auto-correlation between two low-dark counts 
(150 counts/s) single photon counting modules (APD). The S were mounted on a Physics 
Instruments XYZ piezo stage in closed loop operation, (b) Confocal image of 20x20 /im^ 
showing bright spots which corresponds to a color center within a diamond crystal. 



Due to the observed bunching, the center with a ZPL at 749 nm was fit by Eq.(|4]i and 
Eq.(|7]i, which describes three level model. The other three centers, which exhibit two level be- 
havior, were fit by Eq.(l3]i and Eq.dTJi. The deviation from zero of the auto-correlation function 
is attributed to the jitter of the electronics and detectors and residual polarisation dependent 
background. 

The g'iLsiO) of the 744, 749, 756 and 764 nm centers were 0.44, 0.16, 0.2, and 0.09, respec- 
tively. The excited state lifetime of the 744, 756 and 764 nm emitters were 3.8, 3.7 and 13 ns, 
respectively and were obtained by extrapolating the decay rate Ai to zero optical power, refer 
to Eq. |5] see Fig. |5] The fit applied to the second order correlation function for the 749 nm 
emitter contains two exponential decay rates Xi and A2. The extrapolation of these rates to zero 
optical excitation allows the excited state lifetime and shelving decay rates to be determined. 
The decay rate X2 as a function of optical power is shown in the inset of Fig.|5] The resulting 
values of r^^ and are 6.2 MHz and 0.89 MHz, respectively and are much smaller than the 
excited state decay rate, obtained from Ai, r2j=880 MHz (T21 = 1.1 ns). 

The single photon emission count rate as a function of optical power Popt is shown for each 
emitter in Fig|6] The measured count rates (corrected for the background) are given by the sum 
of the counts on the two APDs in the HBT setup. For the 744, 756 and 764 nm emitters, the 
saturation curves were fit according to Eq.(|8]l yielding saturation count rates of 2.1 x 10^, 
3.2x10^ and 1.3x10^ counts/s, respectively. The saturation curve of the three level 749 nm 
emitter was fit using Eq.©, with the values of r2i,ri2, ^31 , r23 obtained from the fit to gmeas(T) 
as a function of pump power and the value 77 obtained as discussed in the following text. The 
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Fig. 3. Typical PL spectra recorded from individual CVD diamond crystals as shown in the 
raster scan of the sample in Fig. 1. The peak emission lines centred at 744 nm (red), 749 
nm (blue), 756 nm (green) and 764 mn (black). Similar emission lines were also found all 
together in one crystal (inset). 

single fitting parameter r]QE was then obtained from a least squares fit to the saturation curve. 

Figure 7 shows a direct measurement of the lifetime recorded from the emitters using the 
pulsed excitation at 20 MHz repetition rate. The single exponential fit to the fluorescence decay 
of each emitter, resulted in measured excited state lifetimes of 4.1 ns, 1.4 ns and 14.2 ns for the 
744 nm, 749 nm and 764 nm centres, respectively. Upon inspection of the fluorescence decay of 
the 744 and 749 nm centres there is a fast component of the fluorescence decay which occurs on 
a time scale less than 0.5 ns, this component is attributed to the background within the emitting 
crystal. The lifetime of the emitter is this case is determined by fitting the fluorescence decay 
for times greater than 0.5 ns. The long fluorescence lifetime of the 764 nm centre mitigates this 
effect and the decay can be described well by the single exponential fit. The measured lifetimes 
are in a good agreement with the lifetimes estimated from the CW measurements for the 744 
nm and 749 nm emitters. The discrepancy with the 764 nm centre can be due to the lack of the 
CW ^'^'(7) measurements at excitation powers well below saturation, which affects the fit to a 
zero excitation power. 

In the inset of Fig. 7 is an example of the second order correlation function measurement 
under pulsed excitation from the 749 nm emitting crystal. The peak at t=0 is the probabil- 
ity of having two photons in the same pulse, indicating in this case the presence of only one 
photon(§'^' (0) = 0.17). The observed deviation from zero is due the background luminescence 
from the diamond crystals, which is not negligible in this case. A pulsed measurement is manda- 
tory for quantum optics application since single photons on demand (triggered single photons) 
are required. To reduce the unwanted background for practical application, temporal filtering 
of our single photon source can be implemented. 

The total collection efficiency, 77, of our setup can been measured directly by exciting a 
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Fig. 4. Background corrected second-order auto-correlation function grneas('f')> measured 
with 154 ps coincidence time bin for 300 s at different optical powers for the (a) 744 nm 
line, (b) 749 nm line the (c) 756 nm line and the (d) 764 nm line. The data of the 749 nm 
line was fit by Eq.l l4l7b . while the data of the 744 nm, 756 nm and 764 nm lines were fit by 
Eq.l l3l7l ). The number to the right of the curves correspond to /"opt/^sat 



two level emitter in the pulsed regime. The total collection efficiency can be described by 
J? = i?opt X i?det X i?e X I7cr, where T7opt corresponds to the optical transmittance, Tjdet is the 
detection efficiency, rjg corresponds to the geometrical efficiency, associated to the matching 
between the collection optics numerical aperture and the source spatial distribution and T7cr 
accounts for the crystal geometry, that contributes to the orientation of the excitation dipole of 
the center with respect to the excitation electric field and the emission dipole with respect to 
the collection optics. These are the limiting parameters that reduce the collection efficiency of 
the source down to typically a few percent. 

In the case of a two-level system the total collection efficiency can be measured directly, as 
each excitation pulse generates one emitted photon (without being trapped in the metastable 
state), hence the absence of detected photons is a directly related to T] . In the condition where 
the laser excitation energy is above saturation, the laser pulses temporal separation is longer 
than the typical detectors dead time (50 ns) and the source lifetime (14.2 ns), the total col- 
lection efficiency can be directly measured from the total count rate and the laser repetition 
rate Tj = (j)/R,-ep- As an example, exciting the 764 nm with a laser repetition rate of 10 MHz, 
77 =1.5% is obtained. Similarly, in the CW regime, rf is given by Tj = ^00/(^25). By calculating 
77 for each two level emitter, an average value of Tj=1.3% is obtained. The discrepancy in the 
value of collection efficiency for various two level emitters may be attributed to a variation 
in rjcr due to the unknown dipole orientation of the emitter within the crystal and/or residual 
polarisation dependent background, that can influence the total saturation count rate (j)oo. The 
accuracy of rjcr can be improved by engineering centers with a known dipole orientation and 




0.6 



Fig. 5. Estimated Aj parameters for the centers at 744 nm (red squares), 749 nm (blue 
diamonds), 756 nm(green triangles) and 764 nm (black circles) versus the optical power. 
At the limit of zero optical power the lifetimes of the centers are respectively (rjj )^'=3.8 
ns, 1.1 ns, 3.7 ns, and 13 ns. The data were fit with Eq.JS]). The behavior of A2 versus 
the optical power of the emitters at 749 nm is also shown inset and fit with Eq.l|6); the 
estimated values r§j =6.2 MHz (161 ns), and r23 =0.89 MHz (1.1 /xs), a = 2.5 mW"', 
j8 =3.1 mW-i. 



lower fluorescent background. A detailed analysis of the various parameters of the collection 
efficiency, which is beyond the scope of the paper, is however important to exploit fully the two- 
level system properties for improving measurement accuracy. After identifying the collection 
efficiency 77=1 .3%, the quantum efficiency 77g£=0.24 of the 749 nm emitter is derived from the 
fit to the saturation curve with Eq.(9). 

A summary of the photo-physics parameters for each single emitter is shown in Table [T] 
The ratio between the average off and on periods of the source is given by the ratio r2ilro,\ 



Table 1. Summary of the photo-physics parameters 
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which equals to 0.14 for 749 nm emitter. This value indicates that the probability of transition 
to the metastable state is moderate, as shown by the slight bunching of the g'^' (t) function (Fig. 
4b). Previous determinations associated to NE8 complex in bulk natural diamond Ila ITOlfTTl . 
and CVD nanodiamond lfT2ll showed ^23/^31 of 2.8, 1.6 and 0.8 respectively, justifying much 
lower saturation rates of NE8 with respect to the 749 nm center. In addition, the center possess 
a shorter lifetime than NE8, which also contributes to the higher emission rates despite the 
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Fig. 6. Measured saturation curves and fit according to Eq.© with estimated parameters 
of optical saturation power, P^^i, and single photon count rate at saturation, (j)co, for the 
centers at 744 nm (red squares) with /'sat=311 A'W and 0c,c=2.1 x 10^ counts/s; at 756 nm 
(green triangles) with Psat=500 and (^oo=3.2xl0^ counts/s; at 764 nm (black circles) 
with -Psat=40 /xW and (j>oo=l3xlO^ counts/s. For the center at 749 nm (blue diamonds) the 
saturation curve was fit by Eq.(|9j with the estimated values ?'2i , ''121 ''31 1 ''23 and 77. From 
the fit a quantum efficiency of 77q£=0.24 is obtained. 

moderately low quantum efficiency. 

In general two level emitters possessed a ZPL with a FWHM of ~ 10 nm, whilst the three 
level emitters exhibited a FWHM of 4 nm. The variation of the lifetime (4 ns to 14 ns) of 
the two level system emitters can be attributed to nanocrystals geometry experiencing local 
modification of electric field and/or differences in the atomic structure. 

The single emitting centres in the range 740-770 nm are believed to arise from chromium 
atoms within the diamond lattice. The charge state of the chromium complex and number of 
atoms involved is unknown. The explanation of what leads to two or three level behavior of a 
particular crystal is still under investigation. One may assume that the strain within the crystal 
may have a significant effect. Indeed, it was shown for NV centers that strain can modify the 
transition from being spin conserving to a A system geometry ||23l. Similarly, the strain can 
modify the electronic transition and enhance/suppress the relaxation to a third metastable state. 

5. Conclusion 

A full optical characterization of the photophysics of novel optical centers attributed to 
chromium related defects in sub-micron CVD grown diamond crystals was performed. Novel 
two-level systems with ZPLs at 744 and 764nm, were identified and characterized. The di- 
rect measurement of the total collection efficiency from the two level systems enabled a more 
accurate determination of the quantum efficiency of the three level single emitter at 749 nm 

(77g£=0.24). 

The spectral width of the two level emitters was ^ 10 nm, whilst the three level emitter 
exhibited a FWHM of ^ 4 nm. This coupled with the short excited state lifetimes (1-14 ns) and 
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Fig. 7. Direct measurement of the lifetime for the centers at 744 nm (red squares), 749 
nm (blue diamonds) and 764 nm (black circles), using a pulsed laser at 20 MHz repetition 
rate with 200 ps pulse width. The data were fit with a single exponential. The deduced 
lifetimes of the centers are respectively 4.1 ns, 1.4 ns and 14.2 ns. Inset: Anti-bunching 
measurement recorded from a single emitter at 749 nm under pulsed laser excitation at 40 
MHz and average power of 70 /iW. 

the high brightness (up to MHz regime), brings these diamond centers much closer to an "ideal 
single photon source on demand". The precise engineering of two-level single photon emitters 
may lead to such sources being implemented as a single photon standard. 
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